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ABSTRACT 


Seismic  analysis  of  RIO  BLANCO,  an  underground 
nuclear  explosion  for  gas  stimulation,  was  performed 
using  LRSM,  LPE,  and  array  recordings.  Although  a 
multiple  detonation,  RIO  BLANCO  was  similar  to  other 
western  U.  S.  shots  and  especially  RULISON,  another 
gas-stimulation  explosion  55  km  away.  Magnitude-yield 
relations,  shear  generation,  and  spectral  content  of 


RIO  BLANCO  were  not  untypical  for  underground  explosions. 


A source  function  which  includes  pP  and  P . . was 

r spall 

inferred  from  the  seismic  data.  An  accurate  location 


was  obtained  for  RIO  BLANCO  using  travel-tire  anomalie 
for  RULISON. 
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INTRODUCTION 


RIO  BLANC3  was  the  third  of  a series  of  gas- 
stimulation  nuclear  explosions  in  the  PLOWSHARE  program. 
Basic  site  information  on  this  shot  is  given  in  Table  I. 
The  feature  of  interest  was  the  multiple  nature  of  the 
shot--actually  :hree  simultaneous  and  closely  spaced 
detonations.  A previous  gas  stimulation  shot,  RULISON, 
was  located  55.9  km  to  the  southeast  of  the  RIO  BLANCO 
site.  RULISON  was  not  a multiple  shot;  therefore, 
comparisons  of  its  signals  with  those  of  RIO  BLANCO 
should  help  to  elucidate  what,  if  any,  effects  the 
multiple  detonation  had  on  seismic  signals. 

In  this  report  we  examine  several  seismic  aspects 
of  the  RIO  BLANCO  shot,  including  location,  magnitude, 
source  function,  and  shear-wave  generation.  Comparisons 
are  made  with  RULISON  especially  and  with  nuclear  explo- 
sions and  earthquakes  in  the  western  United  States  in 
general.  The  report  on  RULISON  by  Lambert  and  Ahner 
(1972)  is  the  basic  reference  for  this  report;  it  will 
often  simply  be  referred  to  as  "the  RULISON  report" 
hereafter. 


-1- 


TABLE  I 


RIO  BLANCO  Detonation  Information 


Date:  17  May  1973 

Time:  16:00:00.1 

Location:  39°47'  34.8"N,  108#21’ 59.6MW 

Surface  Elevation:  6630  ft.  (2010  n) 

Shot  Depths:  5840,  6230,  and  6690  ft. 

(1780,  1900,  and  2040  m) 
(aligned  vertically) 

Uphole  Time:  .527  sec. 

Yield:  30  kt  for  each  of  three 

detonations 

Medium:  Sandstone  (Fort  Union  and 

Mesa  Verde  formations) 
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DATA 

The  data  source  used  for  this  report  on  RIO  BLANCO 
comprised  ten  portable  LRSM  (Long  Range  Seismic  Measure- 
ments) stations,  the  three  on-line  VLLA  arrays  (LASA, 

ALPA,  and  NORSAR)  , the  ten  Long- Period  Lxpenmental  (LPL) 
stations,  and  the  limited  TFO  array.  Magnetic  tape  data, 
cither  analog  or  digital,  were  available  at  the  SliAC  for 
these  sites,  with  the  exception  of  a few  LPL  sites. 

The  ten  LRSM  stations  were  deployed  especially  for 
RIO  BLANCO,  but  usable  long-period  and  short-period  re- 
cordings were  not  available  from  all  sites.  The  LRSM 
recordings  of  RIO  BLANCO  are  displayed  in  Appendix  I. 
Horizontal  components  were  aligned  radially  and  trans- 
versely with  the  shot  location. 

The  arrays  were  all  operational  for  RIO  BLANCO,  but 
short-period  seismometers  at  LASA  and  TFO,  both  within 
8°  of  the  shot,  wc’e  clipped  and  thus  unusable  except  for 
travel-time  determination.  Long-period  signals  were  seen 
at  LASA  and  ALPA  but  not  at  NORSAR.  Ihe  array  recordings 
of  RIO  BLANCO  arc  displayed  in  Appendix  II. 

Data  for  RIO  BLANCO  from  the  ten  LPL  sites  openting 
at  that  time  was  disappointing.  Only  long-period  recor- 
dings are  available  from  these  stations,  and  of  the  ten, 
just  three  stations  detected  surface  waves  (ALQ,  Old),  and 
GIG) . Displays  of  the  RIO  BLANCO  recordings  from  the 
three  stations  are  contained  in  Appendix  III.  Appendix  I\ 
is  a map  showing  the  shot  point  location  and  the  station 
locations  at  which  data  was  recorded  for  RIO  BLANCO. 


The  seismic  data  as  read  from  the  displayed 
recordings  from  RIO  BLANCO  is  presented  in  Table  II. 
Ground  motions  have  been  calculated  using  appropriate 
instrument  response  curves.  Body-wave  magnitudes  have 
been  calculated  using  both  Evernden's  (1967)  and 
Gutenberg  and  Richter's  (1956)  correction  terms;  for 
regional  distances  Evernden's  formulas  for  7.9  km/sec 
and  8.5  km/sec  refractors  were  used  as  indicated,  and 
for  teleseismic  distances  Gutenberg  and  Richter's 
formula  was  applied  as  indicated  by  "G-R".  In  all 
cases  the  same  formulas  as  were  applied  to  RULISON 
stations  (Lambert  and  Ahner,  1972)  were  used  for 
corresponding  RIO  BLANCO  stations.  The  RULISON  report 
discusses  the  choice  of  the  appropriate  formula. 
Surface -wave  magnitudes  have  been  calculated  using 
von  Seggern's  (1970)  distance-correction  term  for 
regional  (A  < 15°)  distances  or  Gutenberg's  (1945) 
distance-correction  term  for  teleseismic  (A  > 15°) 
distances.  Again,  in  all  cases  the  same  formulas  were 
applied  to  RIO  BLANCO  sites  as  were  applied  to  those 
identical  RULISON  sites.  The  body -wave  formulas  are: 

m (7.9)  - 1.21  (log  A/T  + 3.04  log  r)  - 7.55 
(Evernden,  1967) 

rae  (8.5)  « log  A/T  + 2 log  r - 3.27  (Evernden,  1967) 
% (0-  R)  * log  A/T  + B (Gutenberg  and  Richter,  1956) 

where  A is  maximum  zero- to-peak  amplitude  in  milli- 
microns, T is  period,  r is  distance  in  kilometers, 
and  B is  tabulated  by  Gutenberg  and  Richter  (1956). 
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Principal  Phases  for  RIO  BLANCO 











The  surface-wave  formulas  are: 

NT  = log  A/T  + 1.16  log  A + 0.74  A < 15° 

(von  Seggern,  1970) 

M = log  A/T  + 1.66  log  A - 0.18  A > 15° 

(Gutenberg,  1945) 

Where  A now  is  measured  in  millimicrons  peak  to  peak 
on  the  vertical  component.  Note  that  a constant  -0.18 
appears  in  the  last  formula;  this  makes  our  Mg  that 
much  lower  than  would  be  commonly  calculated  by  NOAA 
and  others.  This  constant  appears  when  making  a strict 
conversion  to  vertical  LR  measurement  from  horizontal 
L R measurement  as  originally  proposed  Dy  Gutenberg 
(see  Lambert  and  'diner,  1972). 
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V 


RESULTS 


Location 


Using  the  nearby  RULISON  explosions,  we  ha /e  the 
opportunity  ’;ith  RIO  BLANCO  to  test  the  relative  travel- 
time location  method  in  a new  area.  The  method  is  that 
of  Chiburis  and  Ahner  (1970).  Due  to  the  several  factors, 
mostly  poor  first  arrivals,  it  was  judged  that  only  four 
stations  common  to  RULISON  and  RIO  BLANCO  had  reliable 
enough  data  for  accurate  ? arrival  times.  This  is  the 
minimum  number  required  for  a least-squares  solution  when 
depth  is  constrained  as  was  done  here.  These  times  are 
given  in  Table  III;  the  RULISON  times  were  repicked  by 
the  author  and  do  not  necessarily  correspond  to  those 
the  RULISON  report;  this  was  done  to  insure  utmost  com- 
patibility between  RULISON  and  R.O  BLANCO  signal  start 
points . 


RIO  BLANCO  lies  55.9  km  northwest  of  RULISON.  Using 
the  arrival  times  for  RIO  BLANCO  and  the  anomalies  for 
RULISON  shown  in  Table  III,  the  location  computed  for 
RIO  BLANCO  was  2.04  km  nearly  due  east  of  ground  zero. 

The  final  standard  deviation  of  travel-time  errors  was 
.403  seconds,  and  the  95%  confidence  ellipse  based  on 
the  F-statistic  for  epicenter  location  enclosed  458  km2. 
This  magnitude  of  error  is  remarkably  small  considering 
the  distance  between  the  two  events  and  the  number  of 
stations  used.  The  location  precis  on  is  slightly 
better  than  that  typically  achieved  at  NTS  (Chiburis 
and  Ahner,  1970),  hut  this  may  be  fortuitous  since  the 
sets  of  arrival  times  used  here  are  still  considered  to  be 
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somewhat  uncertain.  However,  this  case  shows  that  the 
minimum  four  stations  may  be  suflicirnt  to  make  a good 
location  when  travel-time  anomalies  are  employed. 

Magnitude 

Taking  all  the  magnitudes,  both  regional  and  tele- 
seismic,  in  Table  II,  RIO  BLANCO  has  an  average 
Ms  * 4.26  and  an  average  mb  = 4.86  with  95$  confidence 
limits  of  +^16  and  ^.28,  respectively.  These  station 
magnitudes  are  listed  in  Table  IV  along  with  RULISON 
magnitudes  wherever  available.  The  RULISON  report 
gives  Mg  = 3.84  ^.15  and  m^  = 4.59  ^.24  for  that  shot 
using  all  available  data  (26  Rayleigh  wave  observations 
and  27  P wave  observations) . A more  accurate  comparison 
of  the  two  shots'  magnitudes  would  use  only  common 
stations  as  seen  in  Table  IV.  With  ju~t  the  nine 
common  stations  for  Ms,  the  surface  wave  magnitude 
difference  between  RIO  BLANCO  and  RULISON  is  0.30;  and 
with  the  eight  common  stations  for  , the  body  wave 
magnitude  difference  is  0.25.  The  yield  of  RULISON  is 
reported  to  be  approximately  40  kt  (Springer  and 
Kinnaman,  1971),  and  that  of  RIO  BLANCO  is  approximately 
90  kt  (Project  Directors'  Completion  Report  D+30  Day:- 
Project  RIO  BLANCO,  Atomic  Energy  Commission).  In 
this  yield  range,  far-ficld  seismic  amplitude  is  both 
predicted  (von  Seggcrn  and  Blandford,  1972)  and 
observed  (Springer  and  Hannon,  1973;  Lvcrndcn,  1970) 
to  be  approximately  proportional  to  yield  for  both 
body  waves  and  surface  waves.  Thus  the  amplitudes 
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TABLE  IV 

Comparison  of  RULISON  and  RIO  BLANCO  Magnitudes 


Station 

RIO  BLANCOS 

RULISON 

h 

RIO  BLANCO0 

RULISi 

BE-FL 

3.88 

3.46 

- 

- 

BP-CL 

- 

4.03 

5.05 

4.58 

CR2NB 

4.56 

4.27 

5.13 

5.02 

HN-ME 

4.58 

4.22 

- 

4.78 

LC-NM 

4.55 

4.11 

4.52 

4.17 

PG2BC 

4.25 

3.89 

4.75 

4.43 

RK-on 

- 

3.48 

4.90 

4.89 

SJ-TX 

4.54 

4.31 

4.30 

4.06 

WH2TfV 

- 

- 

4.67 

- 

WQ-IL 

4.32 

3.99 

- 

4.64 

TFO 

4.26 

4.12 

- 

4.56 

LASA 

4.40 

4.25 

- 

5.06 

NORSAR 

- 

- 

5.56 

- 

ALPA 

3.92 

- 

- 

- 

ALQ 

4.36 

- 

- 

- 

OGC 

3.90 

- 

- 

•• 

CHG 

3.87 

- 

- 

- 

-10- 


r — 


Wk 


for  RIO  BLANCO  should  be  2 h times  those  of  RULISON,  a 
difference  of  0.35  on  the  magnitude  scale.  Our  observed 
differences  of  0.30  and  0.25  are  in  good  agreement  with 
this  prediction,  and  the  lack  of  complete  agreement  can 
most  likel/  be  attributed  to  coupling  and  near-source 
effects,  provided  the  stated  yields  are  accurate.  We 
conclude  that  RIO  BLANCO  behaved  much  as  a single 
concentrated  source  would  insofar  as  magnitudes  are 
concerned.  Considering  the  proximity  of  the  three 
detonations  as  given  in  Table  I,  this  is  not  an 
une/pected  result  for  far- field  observations. 

The  actual  time  delavs  between  the  RIO  BLANCO 
shots  can  be  estimated  as  follows.  The  uphole  time  was 
.527  sec  (Sisemore  and  Toman,  1973);  this  gives  an  uphole 
velocity  of  11,100  ft/sec  (3.38  km/sec)  from  the  topmost 
detonation.  Using  this  velocity  at  shot  level  (though 
a higher  velocity  undoubtedly  applies  there)  the  time 
delays  of  downward  rays  at  the  source  are  .041  sec  and 
.076  sec  for  the  middle  and  topmost  detonation  relative 
to  the  bottommost  one.  Rays  at  any  angle  to  the  verti- 
cal will  be  delayed  less,  by  a factor  of  cos(iQ)  where 
ip  is  the  angle  of  incidence.  For  regional  distances 
this  cosine  factor  should  be  0.7  or  less.  Figure  1 
shows  examples  of  the  effect  on  m^  by  summing  signals 

with  these  time  delays.  The  primary  pulse  is  a synthe- 

_ 7 

tic  one  based  on  the  u>  source  function  discussed  by 
von  Seggern  and  Blandford  (1972).  It  was  created  for 
a 30  kt  source,  a value  T/Q  (travel -time  over  quality 
factor)  of  .54,  and  a LRSM  sho^t  period  instrument. 
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The  two  larger  pulses  are  those  resulting  when  maximum 
delays  are  used  as  stated  above  (appropriate  for  core 
phases)  and  when  0.7  times  these  maximum  delays  are 
used  (appropriate  for  regional  distances)  in  summing 
up  the  three  30  kt  shots  to  simulate  RIO  BLANCO  short 
period  signals.  (The  simulated  signals  had  100  samples 
per  second,  and  so  the  delays  were  rounded  to  the 
nearest  .01  second.)  The  amplitudes  o2  the  summed 
signals  with  delays  are  within  31  of  the  value  that 
would  result  if  the  primary  pulse  were  meioly  multiplied 
by  a factor  of  three  to  represent  a single  90  kt 
detonation.  Also,  no  perceptible  increase  in  period 
occurs.  Thus  the  magnitude  decrease  of  the  three  30  kt 
shots  emplaced  as  for  RIO  BLANCO  relative  to  a single 
90  kt  detonation  should  be  less  than  .02  mfe  unit. 
Howeve.,  this  assumes  no  inelastic  effects  at  the 
source,  whereas  the  RIO  BLANCO  shot  configuration 
probably  led  to  attenuation  of  downgoing  waves  from 
detonations  above  by  the  cavities  and  inelastically 
deformed  surroundings  of  detonations  below.  Our 
conclusions  cannot  be  precise,  but  these  arguments 
suggest  that  the  three  RIO  BLANCO  detonations  should 
not  produce  quite  as  strong  a signal  as  the  combined 
total  of  the  yields  would  predict.  Our  magnitude 
comparisons  of  RIO  BLANCO  with  Rill.  1 SON , considering  the 
error  involved,  docs  not  confirm  this  expectation,  but 
neither  docs  it  refute  this  expectation. 

The  RULISON  report  presented  M^-m^  lor  many 
earthquakes  and  explosions  in  the  United  States,  and 
figure  2 shows  the  RIO  BLANCO  point  (using  all 
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available  data  in  Table  II)  on  that  plot.  The  RULISON 
point  here  represents  all  available  data  in  that  report. 
In  both  eases,  and  also  for  the  other  points  in  this 
plot,  the  appropriate  and  M.  formulas  have  been 
used  for  regional  distances  in  the  western  United 
States.  Although  it  lies  toward  the  upper  extreme  of 
explosions,  RIO  BLANCO  is  apparently  well  separated 
from  wcrtcin  United  States  earthquakes. 

Source  Function 

We  have  already  pointed  out  the  similarity  of  the 
three  RIO  BLANCO  detonations  to  a single  detonation 
insofar  as  magnitudes  avo  concerned.  The  time  delays 
for  seismic  far-ficld  signals  arc  so  small  that  indi- 
cations of  a multiple  source  should  be  non-existent 
in  the  seismic  recordings.  Therefore,  we  do  not  expect 
to  find  multiple-shot  information  by  applying  to.  P wave 
signals  such  techniques  as  spike  filtering,  homomorphic 
filtering,  or  ccpstral  analysis.  However,  the  contri- 
bution of  the  source  function  from  pP  and  spalling 
signals  can  be  studied  with  these  techniques,  and  such 
results  for  RIO  BLANCO  should  be  of  interest  because 
of  its  abnormal  firing  depth  compared  to  most  NTS  shots. 

In  the  RULISON  report,  cepstral  analysis  was 
applied  with  nonuniform  and  indefinite  results  for  the 
LRSM  stations  used;  however,  processing  of  LASA  signals 
indicated  a pP  at  approximately  the  expected  time  and 
also  a later  signal  with  the  same  polarity  as  P which 
may  be  interpreted  as  the  spalling  P wave.  Unfortunately 


A recordings  were  clipped  for  the  RIO  BLANCO  P wave 
ep.tral  analysis  for  the  six  LRSM  stations  with  ade- 
quate S/N  ratios  are  shown  in  Figure  3 This  . 

is  described  by  Cohen  (1969).  Basically  one  exa'min”1  * 
e spectrum  for  nulls  appearing  at  frequencies  appro- 
s or  p or  Pspall  interference  and  the  cepstri- 
ror  a peak  at  the  approximate  time  delays.  For  RIO 
BLANCO.  Sisemore  and  Toman  (1973,  show  the  uphole  time 

1 1 s.l'  ThC-  “ S°  ^ 5h°Uld  te  delaye<i  by  rouRBly 
•1  sec.  Their  accelerograph  data  also  show  the  P 

signal  to  be  delayed  by  .29  sec  at  spal1 

to  DP  anr!  e D f round  zero  relative 

' ° Pspall  should  arrive  roughly  1.4  sec 

after  P on  RIO  BLANCO  recordings.  Using  f„  . n/l 

" 2,"°  f0r  PP  "Ul1  fre<<uenci«,  one  should 

xpect  spectral  nulls  at  0.9,  1.8,  2.7,  etc.  cps  for 

P . and  using  fN  . n/t  („  . i/2,  j/2,  5/2,  etc)  for 

spall  nu  requencies,  one  should  expect  nulls  at 

• . 1.1,  1.8,  2.5,  etc.  cps.  These  nulls  appear  for 

many  of  the  stations  in  Figure  3;  but  the  sum  spectrum 

which  should  smooth  out  individual  path  and  station  ’ 

irregularities  in  the  spectra,  does  not  clearly  contain 

any  of  these  nulls.  All  of  the  individual  stations 

except  LC-NM  show  a spectral  null  i„  the  0. 8-0.9  cps 

range  roughly  correct  for  the  pP  interference,  but 

the  slight  variability  of  the  null  frequency  among 

stations  apparently  leads  to  cancellation  of  this  on 

the  sum  spectrum.  Also  a null  frequency  at  roughly 

‘ ’ tBdlcative  of  P jj,  can  be  seen  on  some 

stations;  but  again  this  is  smoothed  out  in  the  sum 
spectra. 


. 
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Another  technique  available  is  homomorphic 
filtering.  This  is  explained  in  detail  by  Schaeffer 
(1969),  and  it  has  been  employed  with  success  on 
Am’hitkc  Island  test  site  signals  by  Bakun  and  Johnson 
(1973).  Because  the  technique  works  only  on  minimum- 
phase  signals,  we  are  limited  to  data  that  is  beyond 
regional  distances  where  secondary  arrivals  can  make 
the  signal  very  non-minimum  phase.  The  only  P signals 
from  our  data  base  that  were  judged  to  be  amenable  to 
this  processing  were  at  NORSAR  and  WH2YK,  and  a time- 
dependent  weighting  factor  (Schaeffer,  1969)  was 
necessary  to  make  them  even  approximately  minimum  phase. 
The  sampling  rate  for  these  signals  was  20/sec,  and 
the  first  5 sec  of  data  was  extended  to  twice  its 
length  with  zeroes  before  processing.  Figure  4 shows 
the  homomorphic  filtering  results  for  the  two  stations. 
The  weighting  factors  were  ,96^^t  and  .98^0t  for 
NORSAR  and  WH2YK  respectively.  (In  these  figures  the 
complex  cepstrum  is  wrapped  around;  that  is,  the  last 
half  of  the  complex  cepstrum  trace  represents  negative 
time  and  should  be  placed  to  the  left  of  the  start  of 
the  trace.  Significant  amplitudes  at  t<0  indicate 
that  the  signal  analyzed  is  not  completely  minimum 
phase.)  The  complex  cepstrum  was  windowed  with  a box- 
car function  from  -0.8  sec  to  0.8  sec  to  recover  the 
primary  pul*e  and  from  0.8  sec  onward  to  recover  the 
echo  train,  which  hopefully  would  show  spikes  for  pP 

and  P . . with  appropriate  signs  and  delavs  relative 
spall  r 

to  the  primary  signal.  In  both  cases  in  Figure  4,  a 
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negative  excursion  near  1.2  sec  is  in  apparent  rough 
agreement  with  the  1.2  sec  expected  pP  delay.  No 
conclusions  are  possible  for  a PSpan  though. 

A requirement  for  good  results  with  homomorphic 
filtering  is  that  the  primary  signal  and  its  echos  have 
the  same  relative  shape.  But  the  WH2YK  and  NORSAR 
signals  appear  to  have  lower  frequencies  beginning  about 
2 or  3 sec  after  onset  of  the  primary  P.  It  may  well  be 
that  the  spall  signal  is  enriched  in  lower  frequencies 
relative  to  the  explosion  P signal;  the  larger  tem- 
poral and  spatial  scale  of  spalling  relative  to  a 
concentrated  detonation  would  seem  to  require  this 
(Viecelli,  1973).  If  this  is  the  case,  one  would  not 
expect  to  necessarily  see  the  spall  signal  revealed  by 
suppressing  the  short-time  part  of  the  complex  c.epstrum 
to  recover  the  echo  train.  In  fact,  the  change  of 
frequency  content  with  time  on  the  analyzed  record  will 
certainly  distort  the  results  for  the  primary  pulse 
shape  and  pP  delay  time  too. 


We  also  employed  the  inverse  filtering  technique 
described  by  Douglas  et  al.  (1972)  to  NORSAR  and  WH2YK. 
Again,  as  with  homomorphic  filtering,  this  technique 
requires  well- recorded  teleseismic  signals.  The  inverse 
filter  is  that  filter  which  reduces  to  a spike  the 
convolution  of  the  appropriate  instrument  response 
with  an  attenuation  operator  given  by  Futterman  (1962). 

The  Futterman  operator  has  negligible  effect  on  the  phase 
spectrum  of  seismic  signals  at  teleseismic  distances  even. 
T/Q  values  of  .93  for  WH2YK  and  .74  for  NORSAR  were 
used  (since  the  NORSAR  signal  was  beamformed  and  high 
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frequencies  w“re  thus  attenuated  somewhat,  the  value 
of  .74  for  T/Q  is  undoubtedly  higher  than  would  be 
needed  for  a single  instrument  there).  The  inverse 
filtering  results  are  displayed  in  Figure  5.  The 
second  trace  from  the  bottom  is  the  raw  output  when 
the  inverse  filter  is  run  over  the  data  face,  and  the 
bottom  trace  is  merely  the  result  of  using  a low-pass 
filter  on  the  raw  output.  The  bottom  trace  should  be 
an  estimate  of  the  source  displacement- time  function, 
since  instrument  and  path  responses  have  been  removed 
by  the  inverse  filter.  In  both  cases  we  see  the  nega- 
tive excursion  indicative  of  pP  at  approximately  1.2 
sec.  Another  equally  strong  pulse  occurs  positively 
at  approximately  2.0  sec.  According  to  data  at  the 
site  itself,  Pspall  snould  only  be  delayed  by  about 
1.4  sec;  so  either  this  is  not  the  F u signal  or 
that  signal  somehow  had  a delay  approaching  2.0  sec. 
The  only  other  explanation  for  those  strong  positive 
pulses  at  WH2YK  and  NORSAR  other  than  spall  signals 
would  be  secondary  arrivals  along  paths  separate  from, 
but  nearly  identical  to,  that  of  the  primary  P path. 

This  is  a problem  which  our  limited  data  is  unable  to 
resolve . 

Study  of  the  RIO  BLANCO  source  function  has  not 
yielded  precise  results.  We  are  more  confident  of  pP 
evidence  than  Pspall  evidence  since  several  approaches 
have  indicated  a pP  signal  at  approximately  1.1-1. 2 
sec  as  expected.  The  best  evidence  for  P is  in 

the  time  traces  themselves,  where  a lengthening  of  the 
dominant  period,  roughly  2 sec  after  P,  could  reflect 
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Figure  5.  Inverse  filtering  for  RIO  BLANCO  where:  1)  inverse 

filter,  2)  signal,  3)  convolution  of  above,  4)  low- 
pass  filter  output  of  (3),  an  estimate  of  source 
displacement  function. 
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a source  contribution  which  is  more  extensive  in  time 
or  space  than  the  initial  P,  In  addition,  amplitudes 
equal  to  initial  P occurring  3 to  4 ;°c  - "ver  P are 
most  easily  explained  by  a strong  spalling  phenomenon. 


Shear  Generation 


No  direct  S waves  were  identified  on  the  RIO 
BLANCO  recoraings  analyzed  for  this  report.  This 
accords  with  the  fact  that  no  S waves  are  listed  in 
the  data  of  the  RULISON  report.  For  RIO  LLANCO , five 
stations  in  Table  II  have  LDve-wave  amplitudes  reported. 
The  average  LQ/LR  ratio  from  these  five  stations  is 
0.30.  From  the  RULISON  report,  where  nearly  twice  as 
many  stations  were  listed  with  LQ,  one  obtains  a LQ/LR 
ratio  of  0.53.  Thus  the  two  shots  were  within  roughly  a 
factor  of  two  in  their  relative  excitation  of  Love  waves. 
From  all  the  available  LRSM  and  VELA  observatory  data 
for  NTS  and  Amchitka  Island  explosions,  von  Seggern 
(1972)  reports  an  average  LQ/LR  ratio  of  roughly  0.6. 

This  indicates  that  the  Colorado  shots  were  apparently 
somewhat  less  effective  generators  of  Love  waves  than 
the  established  test  sites,  although  there  is  insuffi- 
cient data  to  declare  this  with  any  great  confidence. 


Spectra 


In  Figure  6 wo  show  the  P-wave  spectra  of  RIO 
BLANCO  as  recorded  at  seven  of  the  ten  LRSM  stations 
(no  signals  were  recorded  at  HN-ME  or  BE-FL,  and  WQ-IL 
had  a very  noisy  recording).  These  spectra  are  of  the 
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first  25  sec  of  signal  after  the  P onset,  sampled  at 
20  points  per  second.  Noise  spectra  of  the  25  sec 
preceding  the  P were  used  to  approximately  determine 
the  high  frequency  at  which  S/N  ratio  fell  to  1.0,  and 
the  signal  spectra  were  accordingly  truncated  at  this 
point.  This  was  not  done  at  low  frequencies  since  the 
S/N  ratio  seemed  to  be  above  1.0  to  nearly  the  lowest 
frequencies  sampled.  These  spectra  have  been  corrected 
for  instrument  response  and  are  in  units  of  ground 
displacement.  The  corresponding  RULISON  spectra  were 
computed  for  comparison  in  each  case,  except  for  WH2YK 
where  a recording  was  not  made  for  that  shot.  The 
vertical  scales  of  the  two  shots  have  been  arbitra.ily 
adjusted  for  best  comparison  of  spectral  shape.  The 
spectra  of  the  two  shots  are  very  similar  at  all 
stations  except  LC-NM,  where  low  S/N  ratios  prevailed. 
The  short-period  spectra  :hus  confirm  the  similarity 
of  RULISON  and  RIO  BLANCO  as  seismic  sources. 


Long-period  signals  from  RULISON  and  RIO  BLANCO 
were  near  overlays  in  all  cases,  and  a spectral  com- 
parison was  therefore  not  made  of  these  signals. 
Spectral  ratios  for  LR  calculated  in  the  RULISON  report 
should  apply  to  RIO  BLANCO  within  a few  per  cent. 

That  report  showed  that  RULISON  LR  spectral  content 
was  as  expected  for  explosions. 
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SUMMARY 


Data  from  only  a few  North  American  sites  and  the 
NORSAR  array  were  sufficient  to  locate  and  roughly 
characterize  the  RIO  BLANCO  event.  In  spite  of  the 
multiplicity  of  the  detonation,  RIO  BLANCO  signals  did 
not  differ  in  any  apparent  manner  from  ordinary  explo- 
sions, with  RULISON  as  the  main  comparative  measure. 
Through  homomorphic  filtering,  inverse  ' ltering,  and 
cepstral  analysis,  we  were  able  to  see  the  pP  reflec- 
tion and  possibly  the  spall  impact.  No  direct  shear 
waves  were  identified  for  RIO  BLANCO  and  Love  wave 
generation  was  less  than  that  of  typical  NTS  shots. 
Spectral  content  of  P signals  for  RIO  BLANCO  was  similar 
to  RULISON,  and  LR  signals  were  visually  similar  at 
stations  common  to  both  events. 
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